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NONLINEAR DYNAMIC RESPONSE OF A UNI-DIRECTIONAL MODEL FOR THE 


TILE/PAD SPACE SHUTTLE THERMAL PROTECTION SYSTEM 
By Jerrold M. Housner, Harold H. Edighoffer, and K. C. Park 


SUMMARY 

A uni-directional analysis of the nonlinear dynamic behavior of the space 
shuttle tile/pad thermal protection system is developed and examined for imposed 
sinusoidal and random motions of the shuttle skin and/or applied tile pressure. 
The analysis accounts for the highly nonlinear stiffening hysteresis and viscous 
behavior of the pad which joins the tile to the shuttle skin. Where available, 
experimental data are used to confirm the validity of the analysis. Both ana- 
lytical and experimental studies reveal that the system resonant frequency is 
very high for low amplitude oscillations but decreases rapidly to a minimum 
value with increasing amplitude. Analytical studies indicate that with still 
higher amplitude the resonant frequency increases slowly. The nonlinear pad is 
also responsible for the analytically and experimentally observed distorted 
response wave shapes having high sharp peaks when the system is subject to sinus- 
oidal loads. Furthermore, energy dissipation in the pad is studied analytically 
and it is found that the energy dissipated is sufficiently high to cause rapid 
decay of dynamic transients. Nevertheless, the sharp peaked nonlinear responses 
of the system lead to higher magnification factors than would be expected in 
such a highly damped linear system. 


INTRODUCTION 

The space shuttle orbiter thermal protection system consists of ceramic 
tiles bonded to thin nylon felt pads, known as strain isolator pads, which are 
composed of thousands of intertwined nylon filaments. The pads, in turn, are 
bonded to the aluminum skin (substrate) of the shuttle orbiter. 

During a mission, tile/pad combinations experience dynamic loads arising 
from acoustics, structural vibration, and transonic shock. As a consequence, 
the pad experiences many cycles of loading of varying magnitudes. Experiments 
such as those described in reference 1 have shown that as the pad is cyclically 
loaded and unloaded, hysteresis loops occur in a stress strain characterization 
of the material. Furthermore, these loops creep as a function of stress level 
and number of cycles. The creep of the loops eventually becomes very small with 
each additional cycle, but its effect is to produce a highly nonlinear hardening 
pad material which is quite soft at low stress levels and considerably stiffer 
at higher stress levels. As shown in reference 2, for tiles under static loads, 
the nonlinear pad material properties after cycling significantly affected tile/ 
pad behavior, producing in many cases considerably higher tile/pad through-the- 
thickness interface stresses than before cycling. The nonlinear pad properties 
can also be expected to significantly affect the response of the system under 
dynamic loads. 



For a large class of tile/pad combinations (those for which the tile and 
pad centroids lie along a line normal to the tile surface), a single-degree~of- 
freedom model of the system is sufficient to characterize the response of the 
system under uniform applied dynamic pressure or motion. Such a model can pre- 
dict dynamic tile/pad through-the-thickness interface stresses which are known 
to be critical to the integrity of the system. 

The purpose of this paper is to present the development, verification, and 
application of an analysis for predicting the nonlinear single-degree-of-freedom 
response of the tile/pad/substrate thermal protection system to sinusoidal and 
random applied pressure on the tile and/or imposed motion of the substrate. 

Inasmuch as the mathematical modeling of the pad material is critical 
to the analysis, the modeling will be discussed including the following! 
assumptions; modeling of the observed stress-strain behavior; damping law 
development; and effective percent of critical damping. Verification of the 
analysis is accomplished through comparison with experiment where available. 

In addition, results for the following will be presented: 

(1) Variation of resonant frequency with amplitude of imposed sinusoidal 

substrate acceleration and applied tile pressure 

(2) Variation of magnification factor with amplitude of imposed sinusoidal 

substrate acceleration and applied tile pressure 

(3) Response to random substrate motion. 


MATHEMATICAL MODELING OF THE SPACE SHUTTLE TILE/PAD 
THERMAL PROTECTION SYSTEM 
General Discussion 

As shown in figure 1, the space shuttle tile/pad thermal protection system 
consists of ceramic tiles bonded to the space shuttle through a nylon felt pad. 
In deriving a mathematical model for the system several assumptions are made. : 
For the tile itself, it is assumed, as in reference 2, that 

Assumption (1) : Each tile behaves as a rigid body. 

For the pad, which has a complex nonlinear material behavior simplifying assump- 
tions are made based upon observations of its behavior under cyclic load. 


Conditioning of the Pad Under Slow Cyclic Load 

Material property tests such as those of reference 1 revealed, as displayed 
in figure 2, that under slow cyclic loading of a prescribed amplitude the pad 
exhibits hysteretic behavior. The hysteresis loops thus formed creep as a 
function of load amplitude and number of cycles. However, it is observed that 


2 



the creep of the loop eventually becomes very small with additional load cycles 
of the same or lower amplitude. This process which produces a quasi-stabilized 
material will be referred to as "conditioning." Provided the pad has not 
experienced a load amplitude higher than the one conditioned at, it would appear 
that the creep of the hysteresis loops can be safely neglected during a short 
time analysis. Hence, 

Assumption (2): It is assumed that the pad material has been conditioned at a 

load amplitude which is not exceeded during the period of 
time to be analyzed. 

Analytical responses of the system have indicated that the system is not sensi- 
tive to small additional creeping. As a consequence of Assumption (2), the 
analysis may be carried out over a short period of time (e.g., a few hundred 
cycles) and then restarted using material data corresponding to additional 
cycles at the calculated stress amplitudes which meet or exceed the 
"conditioning" amplitude. 


Pad Through-the-Thickness Stress 

For reasons discussed later in this section, it is assumed that 

Assumption (3): The total through-the-thickness pad stress is the algebraic 

sum of strain-rate independent and strain-rate dependent 
stresses denoted as Oq and O]), respectively, so that the 
total pad stress is 

a = Os + Cd (1) 

Strain-rate independent pad stresses, Og.- It can be observed in figure 2 

that the slow load-rate hysteretic conditioned loop involves Irrecoverable and 
recoverable deformations. The sharply pointed character of the loop indicates 
the presence of static coulomb frictional damping which is strain-rate inde- 
pendent. This is physically reasonable since filaments are sliding over one 
another during loading or unloading. It is, therefore, assumed that 

Assumption (4) : The sIovt load-rate experimental curves account for all 

strain-rate independent stresses, Og, in the pad, both 
elastic and inelastic. 

Figure 3 illustrates typical loading/unloading paths from different strain-rate 
independent states (dg, e) within and on the conditioned loop envelope, where 
e is the through-the-thickness pad strain. 

Based on the observation of pad behavior similar to that exhibited in 
figure 3 and as a consequence of Assumption (2), it is further assumed that 

Assumption (5) : Irrespective of loading or unloading path, the pad strain- 

rate independent state (Og, e) cannot lie outside of the 
conditioned path hysteretic loop envelope. 
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Within the envelope, it is assumed that 

Assumption (6): Loading/unloading paths may be curve fitted approximately by 

appropriate scaling and translation of portions of the 
envelope. 

This assumption is reasonable since the envelope itself is composed of a loading 
and unloading path between the extrema of the conditioning cycle. Details on 
the curve fitting are given in Appendix A. 

Strain-rate dependent pad stresses, O-q.- It is known (see, for example, 

ref. 3), that for a driven linear spring/mass system, the energy dissipation due 
to coulomb strain-rate independent damping is not sufficient to produce bounded 
oscillations at resonant frequencies. Inasmuch as physical oscillations are 
always bounded, it would appear reasonable to postulate that an additional 
energy dissipation mechanism which is strain-rate dependent is present in the 
pad. Indeed, dynamic tests indicate that this is the case. 

It is assumed that the strain-rate dependent stress, 0-qj is given by a 
power law^ namely. 


Assumption (7) : 




(hpe) 


2TT£ 


(q-1 


( 2 ) 


In equation (2), £ is the pad through-the-thickness (z-direction) strain rate, 

m is the tile mass, hp is the pad conditioned thickness, Ap is the tile/pad 
contact area, is a reference frequency, and CjQ are damping parameters to 

be determined from experiment. A single term representation of Oq is used 
rather than a sum of such terms since it is highly desirable to keep the num- 
ber of parameters to be validated from tests to a minimum. The evaluation 
of these parameters is considered in a subsequent section. The existence of 
allows the total stress-strain state (o, £) to lie outside the envelope. 


Integration of the Equation of Motion 
Equilibrium of the tile in the z-direction is satisfied provided 


mw + Ap(ag + ap) = P(t) (3) 

I 

where w is the tile acceleration in the positive z-direction, P is the 
applied tile transient load, 0^ is given by equation (2), and Og is provided 
by tracing out the loading/unloading path which depends upon prior stress-strain 
history where the pad strain is given by 


i 

Examination of such power damping laws in otherwise linear systems is 
contained, for example, in reference 4. 
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e = (w - Wg)/hp 


(4) 


in which w and Wg are the transient displacements of the tile and substrate, 
respectively, in the positive z-direction. 

o Equation (3) is integrated explicitly by using the following recursive 

algorithms at the ith time step: 

Wi+1 = Wi + At|Pi - Ap(0s^ + ODjJ^m (5) 

wi+i = Wi + At Wi+i (6) 


where (•) over any quantity indicates differentiation with respect to time and 
At is an appropriate time step. The initial conditions for equation (3) may 
be specified on e, and £, If the initial state of (a,e) lies outside 

the envelope, then e must initially be nonzero. 


Selection of Strain-Rate Dependent Damping Parameters 

General discussion .- In a linear spring-mass-dashpot system, the response 
amplitude due to sinusoidal oscillations of the base is significantly affected 
by the viscous damping coefficient whereas the resonant frequency is not as 
strongly influenced by the damping coefficient. This is also generally true of 
nonlinear systems except that in addition, the amplitude and resonant frequency 
have a nonlinear dependence on the imposed substrate motion or applied tile 
pressure. In order to evaluate the two damping parameters of equation (2) it 
seems appropriate to select these based upon a comparison of analytical and 
experimental response amplitudes due to imposed sinusoidal motions of the sub- 
strate or applied tile pressure oscillations at the resonant frequency of the 
system. Inasmuch as it is anticipated that the majority of pad stress on the 
space shuttle will be due to substrate motion, the damping parameters are based 
on experimental correlation for imposed sinusoidal substrate motion, where the 
experimental results are presented in references 5 and 6. Once the parameters 
are selected, the system response under sinusoidal applied tile pressure and 
random substrate motion are examined. 

Sinusoidal substrate motion .- Consider the case of P = 0 and an imposed 
sinusoidal substrate motion given by 


Wo = 


gG 


sin 


( 7 ) 


where a is a prescribed amplitude, G is the acceleration of gravity (e.g., 
if Wg is in inches, G is 386 in/sec^), and is the imposed circular 
frequency. 
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The resonant frequency, is taken as that value of ^ which yields 

the peak amplitude response (resonant amplitude) occurring over the frequency 
range of interest (usually 0 to 1000 Hz) . Figure 4 illustrates a typical 
variation with of response amplitude expressed as magnification factor fg, 
where 

fg = w/aG (8) 

For a square LI 900 15.24 x 15.24 x 9.525 cm tile on a 12.7 x 12.7 x Q,406 cm 
pad (6 X 6 X 3,75 in. tile on a 5 x 5 x 0.160 in. pad), figures 5, 6, and 7 , 

present families of analytical curves which show the variations of resonant 
magnification factors with a for different selections of the damping 
parameters C and q. The figures serve to show the effect of each parameter 
on the shape of the curves. Available experimental values are also displayed. 
Within realistic limits it appears that an appropriate selection of the 
parameters is 

C = 0.15, q = 2 (9) 

where the reference frequency, flj., is 100 Hz. 

This selection is guided by the following considerations: 

(1) The results displayed in the figures indicate that irrespective of q 
a value of ? can be found such that the experimental data scatter is reason- 
ably approximated. Hence for q = 1, ? - 0.2; for q = 2, C “ 0.15; for q * 3, 
C ^ 0.06. 

(2) The lowest integer exponential damping law which fits the data scatter 
should be selected; however, at present, experimental data are available only 
for a 5 35. 

(3) In the absence of available experimental data, a criterion for CX > 35 
is required to make a reasonable selection for q. This criterion is taken to 
be that fg must remain bounded as a goes to infinity, in particular 

^^s^max (10) 

where (^s^max chosen to be the maximum value of fg observed among the 
available experimental data. 

(4) The criterion of equation (10) eliminates q = .1, a linear viscous 

damper, as a selection since fg for values of C, which pass through the data 
scatter at a ^ 35, increases above ^(fs)max a increases. 
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(5) The lowest integer exponential damping law which produces a good data 
fit at a ^ 35 and satisfies equation (10) has q and ? as prescribed by 
equation (9). 

(6) Additional experimental data for a > 35 are required to determine 
a better selection. 


COMPUTER CODE 

A computer code which solves equation (1) using the time integration 
algorithm of equations (5) and (6) has been developed. The code traces out the 
loading/unloading curves which are approximated according to the relations pre- 
sented in Appendix A. Relevant response quantities such as pad stress and 
strain, tile and substrate accelerations are output. A short user's guide to 
the code is provided in Appendix B. 


APPLICATIONS 

Sinusoidal Substrate Motion 

Figure 6 provides the variation of the magnification factor, fg, with 
substrate acceleration, a, at the resonant frequency for a LI 900 
15.24 X 15.24 X 9.525 cm tile on a 12.7 x 12.7 x 0.406 cm pad 
(6 X 6 X 3.75 in. tile on a 5 x 5 x 0.160 in. pad). The calculated resonant 
frequency variation with a is given in figure 8 along with corresponding 
experimental results. Agreement is good, with both analysis and experiment 
exhibiting similar trends. 

With increasing substrate amplitude of motion, analysis and experiment show 
a rapid decrease in frequency (effective dynamic softening) followed by a slow 
increase in frequency (effective dynamic hardening). The trend can be under- 
stood by considering the material behavior of figure 3. At low amplitude, the 
material follows a stiff loading/unloading loop (loop A), with little or no 
portion of it on the low modulus region. This indicates that for low level 
oscillations the friction force between fibers accounts for the entire load 
carrying pad mechanism. As the amplitude is increased, a larger portion of the 
loop contains the low modulus region (loop B) , and so the frequency decreases. 
With further increase in amplitude, the frequency starts increasing again as 

portions of the loop begin to include the higher modulus region at raised stress 
levels (loop C). ' 

Examination of the analytically predicted and experimentally observed 
steady state tile/pad interface stress, O, as a function of time reveals the 
highly nonlinear nature of the response. Figure 9 shows <^(t) near the reso- 
nant frequency when a = 30. The agreement between results from analysis and 
test is very good; both show a highly nonlinear response with a sharp stress 
peak. (A linear analysis would give a purely sinusoidal response.) The peaks 
seem to be due to the tile acquiring a high velocity during its travel over the 
soft material range which causes it to overshoot into the stiff material range 
thus producing high stresses which result in the high magnification factors of 
figures 5 to 7. 
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Random Substrate Motion 


Consider again the case of a LI 900 15.24 x 15.24 x 9.525 cm tile on a 
12.7 X 12.7 X 0.406 cm pad (6 x 6 x 3.75 in. tile on a 5 x 5 x 0.160 in. pad), 
where the substrate undergoes random accelerations. Figure 10 presents a time 
slice of the random substrate motion, analytical tile response, and experi- 
mental tile response where the random substrate motion has the power spectral 
density given in figure 11. Correlation between experiment and analysis is 
good in spite of the variability which exists in pad properties from specimen 
to specimen. 

Applied Sinusoidal Tile Pressure 

Using the values of ^ and q established on the basis of imposed 
sinusoidal substrate motion, the system response to sinusoidally varying 
applied tile pressure 


P = + P2^T 


( 11 ) 


is examined for a LI 900 15.24 x 15.24 x 9.525 cm tile on a 
12.7 X 12,7 X 0.406 cm pad (6 x 6 x 3,75 in. tile on a 5 x 5 x 0.160 in. pad) 
where is the tile surface area and 00 is the circular frequency. Fig- 

ures 12 and 13 illustrate the response of the system when Pi = 0. In fig- 
ure 12 the variation of resonant frequency with magnitude of P 2 is presented 
while figure 13 shows the variation of magnification factor, fp, with magni- 
tude of P 2 * In the case of applied tile pressure, fp is defined as 


fp = + P2) 


( 12 ) 


It is evident from a comparison of these figures with figures 5 to 8 that the 
response trends with both sinusoidal substrate and tile pressure excitation 
are similar. 


Application of Analysis To Determine Effective 
Pad Damping Properties 

General discussion .- It has been shown herein that the tile/pad system is 
highly nonlinear. Nevertheless, it may be desirable to establish linear models 
which can provide reasonable results for certain dynamic characteristics of the 
system. To this end, an equivalent linear damper is sought by using the non- 
linear analysis to evaluate an equivalent linear percent of critical damping 
for the pad. 

Five procedures are employed to evaluate the pad percent of critical 
damping. It is shown that different values result depending on the procedure 
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employed. However, for the present need of evaluating pad stresses, a pro- 
cedure based upon stress magnification factor appears to be most suitable. 

Before proceeding with the evaluation, it is useful to review some basic 
linear behavior results associated with base driven, force drive, and free 
oscillations of the system. 

Linear system .- The response of a linear spring/mass/dashpot system with 
spring and dashpot in parallel is well understood. The damping force is pro- 
portional to e and the critical damping coefficient is defined as 


c 


cr 


2x^0 


(13) 


where 0)^ is the resonant frequency without the dashpot. 

Case of an imposed sinusoidal substrate motion: For the case of an imposed 

sinusoidal substrate motion, the magnitude of the mass (tile) motion is given by 


w 


= fa Wc 


(14) 


where near resonance may be approximated by 


fs = ^ n + 


(15) 


in which ^ is the critical damping ratio given by 


e = 


Ccr 2mc0o 


(16) 


The energy dissipated per cycle, Ep, in the case of the sinusoidally 
substrate driven system is given by 


where 


Jd_ 




(17) 


E 


K 



(18) 


is the maximum kinetic energy of the mass during any cycle, and fg is given 
by equation (15). 
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Case of an applied tile pressure: For the case of an applied tile 

pressure, the magnitude of the mass motion is given by 

|w| = fplP|/k (19) 

where k is the linear spring stiffness constant and fp near resonance may 
be approximated by 



( 20 ) 


The energy dissipated per cycle in the case of the sinusoidally force 
driven system is given by 



( 21 ) 


where is given by equation (18) . 

Case of a force free system: For the case of a force free system, in which 

the mass is given an initial displacement and allowed to decay exponentially to 
zero displacement, the critical damping is given by the "log decrement" as 


e = - ^ ln(wt=T/'^t=0) 


( 22 ) 


where 


T = 27T/tOQ (23) 

is the response period. 

Nonlinear system . - For the nonlinear tile/pad system, equivalent values 
of ^ are calculated on the basis of the following five procedures: 

(1) Magnification factor in substrate driven system 

(2) Energy dissipation in substrate driven system 

(3) Magnification factor in tile pressure driven system 

(4) Energy dissipation in tile pressure driven system 

(5) Log decrement 
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where Cj[ denotes a value based on the ith procedure. In a linear system, 
would be independent of force or base acceleration magnitude and be the 
same for all five procedures. However, in the nonlinear pad/ tile system, this 
is not the case. 

Consider the case of a LI 900 15.24 x 15.24 x 9.525 cm tile on a 
12.7 X 12.7 X 0.406 cm pad (6 x 6 x 3.75 in. tile on a 5 x 5 x 0.160 in. pad). 

In figure 14, and ^2 are shown for this case where found from 

equation (15) and the results of figure 6 with C = 0.15 and q = 2 while ^2 
is found from equation (17) where Ej) is calculated by the formula 

% " L ApO(w - Wg)dt (24) 

in which t^ is any time large enough so that the system response is in a 
steady-state condition, w and O are provided by the numerical solution of 
equation (3), and Wg from the differentiation of equation (7). 

The results indicate that except for low substrate accelerations, pro- 
cedures ( 1 ) and ( 2 ) give nearly identical results of about 16 to 18 percent of 
critical damping for moderate to high substrate accelerations. At low accelera- 
tion levels, it would be best to use inasmuch as magnification factors are 

of critical importance to the integrity of the system. 

The results for ^3 and ^ are given in figure 15 where the value of ^3 
is calculated using equation ( 20 ) and figure 12 while is found from equa- 

tion (21) . Unlike the substrate driven system, the pressure driven system 
results in vastly different critical damping ratios between procedures ( 3 ) 
and (4), based on magnification factor and energy dissipation, respectively, 
for all amplitudes of oscillating pressure. Moreover, neither procedure (3) 
nor (4) agrees well with (1) or (2) . Procedure (3) results in damping ratios 
smaller than those of procedure ( 1 ) since the magnification factors under 
applied sinusoidal pressure are larger than those under applied substrate 
motion. On the other hand, procedures (2) and (4) give comparable values for 
peak critical damping ratios, but procedure ( 4 ) does not result in the drop off 
of the damping ratio seen in procedure ( 2 ) for increasing input energy. 

The results of ^3 are given in figure 16 where the analysis has been used 
in conjunction with equation (22) . Procedure (5) results in a monotonically 
increasing critical damping ratio with increasing initial strain and thus yields 
values seen in all the other four procedures. Use of equation (22) implies a 
full cycle determination of log decrement whereas reference 4 suggests the use 
of a half cycle determination when q > 1. However, since the pad properties 
lead to generally larger responses on the tensile half cycle than the compres- 
sive half cycle even when there is no viscous damping, a half cycle determina- 
tion would produce misleading results whereas a full cycle determination would 
not . 

Two conclusions are drawn from the comparison of these five procedures. 
First, the system is too nonlinear to permit an equivalent linearized damping 
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evaluation that allows for all possible situations. The nonlinear behavior 
gives rise to high stress peaks and thus high magnification factors even in the 
presence of considerable energy losses due to pad dissipation mechanisms. 

Second, if a linearized model is desired, the critical damping ratio should be 
chosen so as to correlate with the nonlinear results which are of interest to 
the analyst. In the case of the space shuttle, this implies using about 16 per-- 
cent or 10 percent critical damping for imposed substrate motion or applied tile 
pressure, respectively. These numbers are based on magnification factors, pro- 
cedures (1) and (3), respectively. 

Energy dissipation due to friction alone .- By taking the damping parameter, 
C, equal to zero, the energy dissipated per cycle as a consequence of strain- 
rate independent frictional mechanism in the pad may be examined. This is shown 
in figure 17 for imposed sinusoidal substrate acceleration. Comparison with the 
results of procedure (2) in figure 14 indicates that the energy loss due to 
friction alone is about 25 to 30 percent of the total energy loss with the 
remainder of the loss due to nonlinear viscous damping. 


CONCLUDING REMARKS 

A uni-directional nonlinear dynamic analysis and associated computer code 
for the space shuttle tile/pad thermal protection system have been developed and 
examined for imposed sinusoidal and random motions of the substrate and/or 
applied tile pressure. The analysis accounts for the nonlinear stiffening, 
hysteresis, and viscous behavior of the pad material and is in reasonably good 
agreement with experiment for resonant frequency and response magnification 
under imposed sinusoidal substrate motion as well as imposed random substrate 
motion. Application of the analysis reveals a highly nonlinearly behaving 
system whose resonant frequency is very high (over 500 Hz) for low amplitude 
oscillations, but rapidly decreases to a minimum value with increasing ampli- 
tude from which it increases slowly with further increase in amplitude. The 
nonlinear pad also causes distorted response wave shapes having high stress 
peaks when the system is subject to sinusoidal loads. Furthermore, energy 
dissipation in the pad per cycle leads to relatively high ”log decrement" 
damping and thus rapid decay of dynamic transients. Nevertheless, the sharp 
peaked nonlinear responses lead to higher magnification factors based on peak 
responses than would be expected in an equivalent highly damped linear system. 
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APPENDIX A 


CURVE FIT OF PAD MATERIAL PROPERTIES 
Conditioned Envelope 

As a consequence of assumptions (5) and (6) in the main body of the paper, 
a curve fit of the conditioned hysteretic loop is first sought. To this end 
the envelope is divided into six segments labeled sequentially in figure 18. 
Segments (1) and (2) are straight lines whose equations are defined in terms of 
measured values indicated in figure 18 as 


si " ®b + 

(Al) 

s2 ~ ^a ^z^ 

(A2) 


Parts (3) to (6) are defined by fifth order polynomials as 


^si ^iO ^ ^il^^ ^il^ ^13^^ ^11^ 

+ Ci5(e - i = 3, 4, 5, 6 


(A3) 


where the four coefficients for the ith curve are chosen such that the ith curve 
passes through the measured points of , £ij), j = 1,2,3, on part (1) of the 

envelope and has a slope equal to E^ at (a^i, £ii) . Since parts (4) and (5), 
and (3) and (6) join at their respective end points. 


*^43 " '^535 ^43 - ^535 <^33 “ ^635 ^33 " ^63 (^^) 

and Eg is given by 

" (^51 " ^31^/(^51 “ ^3l) (A5) 

while 0 ^ and 0^ are given by 

^a ^61 " ^z^61 

(A6) 

'^b " ^31 ~ ®z^31 
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and 


'^41 ~ ^61 ^ ^51 ~ ^31 
^41 " ^61 ^51 ■ ^31 


(A7) 


Consequently, nineteen independently measured points are required to charac- 
terize the envelope curve fit used herein. 


The coefficients of equation (A3) may now be determined as 


*^i0 ■ '=^il 


Cil = E, 


(A8) 


and 


r > 


'i3 


- 


“"iS 


(ii2 - £ii) (e^, - e,-i) 


i2 "-il^ 


^^i3 ~ ^il^^ (^i3 " ^il^^ 


-1 


^i2 ‘^il ^z^±l 


^i3 ■ ^il " 


(A9) 


Paths Within the Envelope 

General discussion .- With the envelope defined by equations (Al) to (A3), 
the constraint associated with assumption (5) can be met. Paths within the 
envelope by assumption (6) follow trajectories similar to the envelope itself. 
These trajectories are designated as positive or negative paths depending on 
whether the strain rate, e, is positive or negative along the path. For 
example, on the envelope, paths (1), (3), and (5) are positive whereas (2), (4), 
and (6) are negative. Paths within the envelope are determined by the state of 
stress-strain existing when e changes sign. The approximations to the posi- 
tive and negative paths are now described. 


Positive paths .— Positive paths within the envelope are generally approxi- 
mated by a curved_portion, or a straight portion, depending on whether the 
reversal strain when e last changed sign from negative to positive is 

less than or greater than the upper cutoff strain, iy. However, when a posi- 
tive interior path intersects a positive envelope path, the envelope path is 
followed . 

Reversal strain less than upper cutof f_strain: When the reversal strain, 

£a» is less than the upper cutoff strain, e^, the positive path is approximated 
by a fifth order polynomial similar to segment (3) but passing through three 
points designated as A, B, and C and with slope at point C. The state of 

stress-strain at these points is as follows: 
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Point A has the state of stress-strain (5^, e^) at which e changes sign 
from negative to positive. 

Point B has the state of stress-strain (Og, Eg) proportioned from the 
envelope curve (3) and given by 



^32^z 


■ ^32 


where 


"Z ^ ^ ^ , 

^B ■■ 


^Bp - ^Cp 


^33 


- Or 


(a. - Or*) 


^C* ~ ^31 ^^A “ ^31^®z 

^Bp ^32 ^32^z 

^Cp ° ^31 

^31 " ^(£33 - £33) 

^32 = ^(£33 - £32) 


(AlO) 


(All) 


Point C has the state of stress-strain (Oq, Oq) given by 


= *^31 + (^A ~ ^31) ^z ~ ^31’^z 

(A12) 

= ^A " ^31 


The factor Y ~ 1 is introduced because the paths at low stress levels cannot 
be accurately determined from test due to the very soft nature of a low stressed 
pad. Its value is 

Y = ; £61 - £a - £51 (A13) 
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where 


r 




and unity otherwise. Generally y ”>ay be taken as 1/2. 


The coefficients of the polynomial which must pass through points A, B, 
and C can be readily determined from equations (A9) . This is_accoraplished by 
replacing 0 ^;^, 0 ^ 2 > ^il ^A» ^C> respectively, and £^ 3 , £i 2 > ^il by 

^A’ ^B’ ^C* respectively. The form of equation (A3) provides the strain-rate 
independent stress component, Og, on the appropriate path. 

Reversal strain greater than upper cutoff strain: When the reversal 

strain, e^, is greater than the upper cutoff strain, £u> bhe positive path is 
approximated by a straight line segment of slope El equal to the secant slope 
of the previous negative path times a factor 3 which is chosen to agree with 
experimentally observed behavior. 


Negative paths .- Negative paths are similar to the positive paths but 
employ a reversal strain, » when £ last changed sign from positive to 
negative and a lower cutoff strain, £l» When a negative interior path inter- 
sects a negative envelope path, the envelope is followed. 

Reversal strain less than lower cutoff strain: When the reversal strain, 

E^i , is less than the lower cutoff strain, £l, the negative path is approxi- 
mated by a straight line of slope Ey equal to the secant slope of the previous 
positive path. 


Reversal strain greater than lower cutoff strain: When the reversal 

strain, , is greater than the lower cutoff strain, the negative path 

is approximated by a fifth order polynomial similar to segment (4) but passing 
through three points designated as A* , B', C’ and with slope E^ at point C'. 
The relationships providing the state of stress-strain at these points may be 
found from equations (AlO) through (A12) by replacing the subscripts A, B, C, 
Og, and 3 by A' , B', C’, o^, and 4, respectively. The coefficients of equa- 
tions (A9) can then be evaluated from the resulting relationships, and then the 
form of equation (A3) provides the strain-rate independent stress components, 
Og, on the appropriate path. 

Cutoff strains, £y and El'" These strain values determine the regions 

of the enclosed envelope where straight line approximations are employed for 
positive and negative paths, respectively. Inasmuch as considerable computa- 
tional time may be involved in creating polynomial path approximations each 
time the strain rate changes sign it can be advantageous to afford some loss in 
accuracy by using straight line approximations everywhere. By setting 
£u S £33 and £l ^ P53 only straight line approximations are used. On the 


) 
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other hand, by setting Sjj = = 0 curved negative paths with straight line 

positive paths are used when > 0 and curved positive paths with straight 

line negative paths are used when 

Additional restrictions: The tangent slope on any path may not exceed the 

maximum tangent slope on the envelope and must always be positive. In order to 
avoid poor curve fit approximations, should be less than £5^^ and e-^ 

should be greater than £51* 


17 



APPENDIX B 


USER INSTRUCTIONS FOR DYNOTA (DYnamic ^nlinear Tile Analysis) 

General Discussion 

This appendix describes the user input and output of DYNOTA which is a 
computer program based on the uni-directional nonlinear dynamic tile/pad 
analysis in the main body of this paper. The program computes the response of 
the tile/pad system to sinusoidal and random imposed substrate motion and/or 
applied tile pressure. 


Input Description 

The program requires as input the number of cases to be executed, program 
control parameters, tile and pad geometry data, nonlinear viscous damping param- 
eters, loading data, initial conditions, time step, and pad loading/unloading 
data. Following the first record which prescribes the number of cases to be 
executed, each case requires a minimum of 5 records and a maximum of 13 records 
as shown in figure 19. 

The minimum 5 records are sufficient when the user wishes to take advantage 
of the standard loadirig/unloading pad property data built into the code and 
imposed substrate motions and applied tile pressures are sinusoidal in time. 

When imposed substrate motions and/or applied tile pressures are random one or 
two additional cards are required. If the standard pad data are not desired, 
the code allows for limited flexibility in adjusting these data by the use of a 
single additional record or complete flexibility by accepting user supplied pad 
property data on 5 additional records. A description of all records follows. 

Record A1 - This record prescribes the number of cases to be executed as 

NCASE 

according to format 15. 

Record B1 - This record prescribes control parameters 

ISUB, IPRESS, MAXSTP, IRATE, IPRINT, IPLOT 

according to format 615 where 

ISUB = 0, no substrate motion 
= 1, substrate motion 

IPRESS = 0, no applied tile pressure 
= 1, applied tile pressure 

MAXSTP = total number of time steps 
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IRATE = 0, selects standard loading/unloading curve fit data for 
pad material properties without adjustment 

= 1, expands or contracts the standard curve fit data parallel 
to the strain axis by a factor STRPER (defined on the 
D1 record) and moves the lower envelope of the standard 
data up by an amount DSTRS (defined on the D1 record) 

== 2, user supplies loading/unloading curve fit data on 
records D2 through D6 

IPRINT: Controls the values of time at which print output of pad 

strain, pad stress, rms stress, tile acceleration, 
substrate acceleration, and applied pressure are displayed 

= 0, prints only at those times when pad stress is a 
relative maxima or minima 

= 1, prints every time step 

= N, prints every Nth time step 

IPLOT: Controls output plots 

= 0, no plots 

= 1, generates the following plots for random response: 

(a) Imposed random substrate acceleration versus time 

(b) Applied random tile pressure versus time 

(c) Tile response acceleration versus time 

= 2, same as IPLOT = 1 but for pure sinusoidal imposed 
substrate motion and applied tile pressure 

=3, plots out the stress-strain curves being traced during 
entire period of case execution 

Record B2 - This record prescribes file and pad properties as 

AREAP, AREAT, THICK, SM, GAMMA, BETA, CUTOFFl, CUT0FF2 

according to format 8E10.4 where 

AREAP = contact area between pad and tile 

AREAT = surface area of tile on which pressure acts 

THICK = pad thickness after conditioning 
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N 


SM = tile mass 

GAMMA = low stress loading/unloading factor, see equation (All) 
(Default: 1/2) 

BETA = multiplying factor for straight line loading /unloading 

interior paths as described in Appendix A (Default: 2) 

CUTOFFl = iy, see Appendix A (Default: £ 51 ) 

CUT0FF2 = see Appendix A (Default: £ 51 ) 

Record B3 - This record prescribes pad damping properties as 

DAMP, QM 

according to format 2E12.5 where 
DAMP = C (Default: 0.15) 

QM = q (Default: 2) 

Record B4 - This record prescribes applied substrate motion and tile pressure 
AMPS, FREQS, AMPPl, AMPP2, FREQP, GRAV 
according to format 6E10.4 where 

AMPS = a = dimensionless amplitude of sinusoidal substrate accelera- 
tion in multiples of G 

FREQS = S^/2tt = circular frequency of sinusoidal substrate accelera*" 
tion (Hz) 

AMPPl = Pj^ = amplitude of nonoscillatory component of applied 
pressure 

AMPP2 = P 2 = amplitude of sinusoidal component of applied tile 
pressure 

FREQP = u)/2Tr = circular frequency of sinusoidal component of 
applied tile pressure 

GRAV = G = gravitational constant (Default: G = 386 in/sec^) 

Note for Random Response : 

FREQS and/or FREQP set to zero implies Imposed random substrate accelera- 
tion and/or applied random tile pressure, respectively, read from binary 
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tapes, TAPEl and/or TAPE2, respectively. Each record contains one 
data entry. 

Record B5 - This record prescribes initial conditions and time step as 

SIG, EPS, EPSD, DT 
according to format 4E10.4 where 
SIG = initial pad stress 
EPS = initial pad strain 
EPSD = initial tile strain rate 
DT = time step, At 

Record Cl - This record is used only when ISUB = 1 and FREQS =0 to 
prescribe specifications on imposed substrate acceleration 
histories on TAPEl as 

NRl, NSl, TINCl, SCALEl 
according to format 215, 2E10.4 where 
NRl = total number of records on TAPEl 
NSl = number of records at start of TAPEl to be skipped 
TINCl = uniform time increment between records on TAPEl 
SCALEl = scale factor to be applied to TAPEl data 

Record C2 - This record is used only when IPRESS = 1 and FREQP = 0 to 
prescribe specifications on applied tile pressure histories on 
TAPE2 as 

NR2, NS2, TINC2, SCALE2 
according to format 215, 2E10.4 where 
NR2 = total number of records on TAPE2 
NS2 = number of records at start of TAPE2 to be skipped 
TINC2 = uniform time increment between records of TAPE 2 
SCALE2 = scale factor to be applied to TAPE2 data 
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Record D1 - This record is used only when IRATE =1 on the B1 record and 
prescribes adjustments to the standard loading/unloading pad 
property data as 

STRPER, DSTRS 

according to format 2E10.4. See description under IRATE *» 1 on 
B1 record for definitions. 

Records D2-D6 - These five records are used only when IRATE =2 on the B1 
record to prescribe user supplied loading/unloading pad 
property data according to the formulation of Appendix A 
(illustrated in figure 19). These records have the form 


*^a» ^b» 

z 


> 

record 

D2 

^31» ^32’ 

^33’ 

^32’ 

^33’ 

record 

D3 

^41’ ^42’ 

^43» 

^42» 

^43’ 

record 

D4 

£51* £52. 

052 


9 

record 

D5 

^61’ ^62» 

^62 


9 

record 

D6 

each according 

to format 

5E10. 

.4 



Output Description 

Program displays at user selected times pad through- the-thickness stress 
and strain, tile acceleration, substrate acceleration, and applied pressure. 
Also displayed are rms stress, tile acceleration, and imposed substrate 
acceleration. In addition, the program assembles necessary arrays of data for 
post-processing. The user may use these arrays in conjunction with post- 
processing routines such as those which compute power spectral density. Such 
routines are generally available in the program library of most computer 
centers. 
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Figure 3.- Loadii^/unloading paths within a conditioned hysteretic envelope 

[pad thickness = 0.406 cm (.106 in.)] 



Magnification factor. 
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Imposed sinusoidal substrate frequency, Hz 


Figure 4.- Typical variation of magnification factor with imposed sinusoidal substrate 
frequency for a prescribed substrate acceleration amplitude, a = 10, (q = 2, C = -15). ^ 
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Figure 5.- Effect of nonlinear viscous damping parameter, 
on magnification factor when q = 1., 0^=100 Hz. 
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Figure 8.- Variation of tile resonant frequency with imposed 

acceleration amplitude, (q = 2, ? = .15). 
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Figure 9.- Nonlinear wave shape due to sinusoidal substrate acceleration, a - 30. 
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Figure 14.- Effective percent of critical damping on the basis of magnification factor 
and energy dissip^ion/ cycle under imposed sinusoidal substrate acceleration. 
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Figure 15.- Effective percent of critical damping on the basis of 
magnification factor and energy dissipation/cycle under applied 
sinusoidal tile pressure. 
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Figure 16.- Effective percent of critical damping on the 
basis of log decrement, procedure 5. 







Figure 18.- Curve fit for strain-rate independent loading/unloading path envelope. 
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Figure 19.- Schematic of input data deck. 
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